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A method for calculating glass sheet deformation based on a general model at temperatures above the vitrifi¬ 
cation temperature is proposed and allows analytical determination of design parameters of molding equip¬ 
ment elements and their operating conditions. 


Relaxation of deformations and stresses is extremely sig¬ 
nificant for glass rheology. In particular, S. M. Rekhson and 
V. A. Ginzbutg in their studies [1 ] discuss a model describing 
the relaxation of the glass structure, stresses and deforma¬ 
tions (IKhS model) which clarifies the character of modifica¬ 
tions of the parameters mentioned above within the vitrifica¬ 
tion interval. The model uses an arbitrarily set history of 
variation of stress a and deformation 6 in time x. According 
to [1], the total deformation is described by the following ex¬ 
pression 

£(t) = e d (t) + e n (r) + e(0), (1) 

where ejx) is the delayed elastic component of deformation; 
e (x) is the viscous component of deformation; e(0) is the in¬ 
stantaneous (elastic) deformation component. 

A graphic interpretation of these types of deformation is 
presented in Fig. 1. 

According to the theory of solids, delayed elastic defor¬ 
mation is similar to creep. 

In calculation of deformation based on Eq. (1), the most 
complicated is determination of e/x), while the other com¬ 
ponents can be easily calculated from the formulas 

e*(x)=ct(0)y; 

n 


In the general case, according to the data in [ 1 ], the creep 
component of deformation is found from the expression 

E d M = e d (°°)( 1 " A'W)’ ( 2 ) 

where e^oo) is magnitude of the function e^x) at x -»• oo; 
A/.(x) is the deformation relaxation function. 

The magnitude of the function M* is calculated from the 
following formula derived from the recommendations in 
[1-3]. 

A/ E (x) = exp(-(x/x r ) fc ), 

where x r is the relaxation time; b is a constant (b = 0.5) 

In Eq. (2), only parameter e^oo) is unknown. The method 
for determination of this component is based on the geomet¬ 
ric decomposition of function e(x) into its components (Fig. 2). 



%(<») 


e(0) = 


o(0) 

£ 


where o(0) is the initial stress under loading; x is the loading 
or relaxation time; q is the glass viscosity £ is the modulus 
of elasticity in tension. 
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Fig. 3. Temperature dependence of x c . 




Fig. 4. Design diagram for determining ( 0 (g) under loading with a 
distributed load. 


According to Eq. (2), function e/x) approaches a certain 
value £/!)„„„ which is nothing else but e/®). Thus, it can be 
postulated that as soon as e/x) reaches a certain value, func¬ 
tion e(x) depends only on the viscous component e n (x), and 
the contribution of e/x) to the total deformation is limited by 
a certain constant value. The characteristic point in this case 
is e(x c ) at which, at time x c , the values of components e/x) 
and e/x) become equal. 

Assuming that for constant temperature and loading con¬ 
ditions e/®) is a constant value, x c can be found, provided 
that the value of the function 

/(x) = l-M/x) 

for this point is known. This value is determined experimen¬ 
tally (Fig. 3). Then, denoting the value of /(x) as cp we can 
write: 


Fig. 5. Variation of sample deflection depending on tempera¬ 
ture: 1, 2, 3, and 4 ) 620, 630, 650, and 660°C, respectively. 


Using the value obtained, e/®) can be determined: 


( 00 )= 


g (°) T c 

tK1-A/ e (t c ))' 


Since all components of function e/x) are found, total 
deformation e(t) can be determined as well. 

Using the expression obtained for total deformation e(x), 
one can predict the behavior of a glass sample both in defor¬ 
mation and in relaxation of deformation, which is the most 
important moment of sheet glass deformation. 

Consider the application of the given model to bending 
deformation. The loading diagram is shown in Fig. 4. For 
this loading scheme, the equation of the bent beam axis can 
be described by equation [4] 


ifl(Z)=- 



24 E 2 J 



where q is the distributed load; L is the beam length; E 2 is 
the shear modulus; J is the moment of inertia of the beam 
lateral section; Z is a longitudinal coordinate. 

Then, using the expression 

j M (Z ) max 

P dr 2 E 2 J ’ 

h 1 

max 2 p ’ 


x c = x r (- ln(! — cp)) 1 '* - 
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where p is the radius of curvature of the beam middle layer; 
M (Z ) is the maximum torque; h is the beam thickness, 
one can obtain the expression for to(e): 


Fig. 5 shows the variation in the deflection g>(t) of the 
sample depending on the temperature at the loading site. Re¬ 
laxation curves are constructed and calculated with the same 
formulas taking into account relaxation of stress in the hold¬ 
ing time segment. These curves were obtained on the basis of 
the model considered and the sample loading diagram is 
shown in Fig. 4. The calculations were carried out for the fol¬ 
lowing loading procedure: loading of the sample, holding 
time (relaxation of stress), removal of the load (relaxation of 
deformations). The following initial data were used: defor¬ 
mation time of 1 sec; beam length of 0.1 m; beam thickness 
of 0.004 m. 

On the plots in Fig. 5, the elastic component is absent. 
That was done for more accurate rendering of the process. 


This particular component can be determined by the respec¬ 
tive formula. 

Based on the curves obtained, the time and temperature 
parameters of the glass sample deformation can be deter¬ 
mined for prescribed limiting values of residual stress defor¬ 
mation relaxation. 

The method considered will become the basis for deter¬ 
mining the design parameters of molding equipment ele¬ 
ments and their operating conditions. 
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